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ABSTRACT 

o . 

£ — , Context. Several studies of nearby active galaxies indicate significantly higher HCN-to-CO intensity ratios in AGN (e.g., NGC 1068) 

than in starburst (e.g., M82) environments. HCN enhancement can be caused by many different effects, such as higher gas densities 
and/or temperatures, UV/X-ray radiation, and non-collisional excitation. As active galaxies often exhibit intense circumnuclear star 
formation, high angular resolution/high sensitivity observations are of paramount importance to disentangling the influence of star 
formation from that of nuclear activity on the chemistry of the surrounding molecular gas. The tight relation of HCN enhancement 
and nuclear activity may qualify HCN as an ideal tracer of molecular gas close to the AGN, providing complementary and additional 
5_l ' information to that gained via CO. 

Aims. NGC 6951 houses nuclear and starburst activity, making it an ideal testbed in which to study the effects of different excita- 
tion conditions on the molecular gas. Previous lower angular resolution/sensitivity observations of HCN(l-O) carried out with the 
Nobeyama Millimeter array by Kohno et al. (1999a) led to the detection of the starburst ring but no central emission has been found. 
Our aim was to search for nuclear HCN emission and, if successful, for differences of the gas properties of the starburst ring and the 
k^l nucleus. 

Methods. We used the new A, B, C and D configurations of the IRAM PdBI array to observe HCN(l-O) in NGC 6951 at high angular 
,— - . resolution (1"=96 pc) and sensitivity. 

Results. We detect very compact (<50pc) HCN emission in the nucleus of NGC 6951, supporting previous hints of nuclear gas struc- 
ture. Our observations also reveal HCN emission in the starburst ring and resolve it into several peaks, leading to a higher coincidence 
between the HCN and CO distributions than previously reported by Kohno et al. (1999a). 

Conclusions. We find a significantly higher HCN-to-CO intensity ratio (>0.4) in the nucleus than in the starburst ring (0.02-0.05). As 
for NGC 1068, this might result from a higher HCN abundance in the centre due to an X-ray dominated gas chemistry, but a higher 
gas density/temperature or additional non-collisional excitation of HCN cannot be entirely ruled out, based on these observations. The 
compact HCN emission is associated with rotating gas in a circumnuclear disk/torus. 

Key words. Galaxies: individual (NGC 6951) - Galaxies: active - Galaxies: nuclei - Galaxies: Seyfert - Galaxies: starburst 



1. Introduction 

Over the past decade, several studies (e.g., Tacconi et al. 1994, 
Sternberg et al. 1994, Kohno et al. 1999a,b, 2000, 2001, 2005) 
have shown that the HCN-to-CO intensity ratio (=Rhcn/co) 
can be significantly higher in active galactic nuclei (AGN; e.g., 
NGC 1068) than in starburst (SB; e.g., M82) or quiescent re- 
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gions. While surprisingly large Rhcn/co (up to 1 or more) have 
been observed in AGN (e.g., NGC 1068, NGC 1097, NGC 5194; 
Kohno et al. 1999a, Tacconi et al. 1994, Usero et al. 2004), 
much smaller Rhcn/co (<0.3) are detected in pure SB or com- 
posite (AGN+SB) galaxies (e.g., Arp220, M82, NGC 6951; 
Gao & Solomon 2004a,b, Nguyen-Q-Rieu et al. 1992). Inactive 
galaxies have even lower ratios of Rhcn/co<0-1. Many differ- 
ent effects can contribute to increased Rhcn/co in active envi- 
ronments including higher gas opacities/densities and/or tem- 
peratures, non-standard molecular abundances caused by strong 
UV/X-ray radiation fields or additional non-collisional excita- 
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Fig. 1. Integrated HCN(l-O) emission (black contours) overlaid on CO(2-l) (color scale; Schinnerer et al., in prep., Garcfa-Burillo 
et al. 2005) in natural (left, la) and uniform weighting (right;Vo); the CO and HCN emission have been both integrated from 
-200km s _1 to +200km s _1 . We used a uv-taper for CO to match the angular resolution of our HCN data which is by a factor of ~2 
lower, and obtain identical beamsizes. Black contours run from 3<x to Tier (right: 7<x) in steps of l<x=0.06 Jy beam -1 km s (right: 
lcr=0.10 Jy beam -1 km s _1 ); CO(2-l): 3cr=1.0 Jybeam -1 km sT l (right: 3cr=1.3 JybeairT 1 km s _1 ). The black line (left) indicates 
the major axis of the bar (PA=100°) and the grey line the major axis of the galaxy (PA=130°). The black line (right) represents the 
most extreme central velocity gradients (PA=(160+20)°). The (0,0) position is at aj2ooo=20 h 37 m 14.123 s and (5 J2 ooo=66°06'20.09" 
which is slightly different from the phase centre of the observations. 



tion such as IR pumping through UV/X-ray heated dust. Gao 
& Solomon (2004a) have ruled out the latter scenario for large 
scale HCN emission and the lack of any clear correlation be- 
tween the hard X-ray and MIR luminosity in AGN (Lutz et al. 
2004) reduces the significances of IR pumping also at small 
scales. However, Usero et al. (2004) present strong evidence in 
the case of NGC 1068 that the nuclear gas chemistry is dom- 
inated by X-ray radiation from the AGN yielding significantly 
different molecular abundances than in SB or quiescent environ- 
ments (Lepp & Dalgarno 1996; Maloney et al. 1996). Recent 
IRAM 30m observations of several HCN transitions in a sample 
of 12 nearby active galaxies also seem to support a significantly 
higher HCN abundance in AGN than in SB or inactive environ- 
ments, rather than a pure density/temperature or non-collisional 
excitation effect (Krips et al., in prep.). If true, this has a se- 
vere impact on the interpretation of Rhcn/co as a measure of 
the dense to total molecular gas mass fraction in active galaxies 
(e.g., Gao & Solomon et al. 2004a,b) as discussed in Gracia- 
Carpio et al. (2006). The study of nearby active galaxies also 
reveals the limitations of Rhcn/co as a unique diagnostic in dis- 
tant sources whose starburst and AGN components cannot be 
separated. 

NGC 6951 is an active galaxy of Hubble type SAB(rs)bc at 
a distance of 24 Mpc (Tully 1988); its active nucleus is classi- 
fied as a transition object between a LINER and a type 2 Seyfert 
(Perez et al. 2000). In addition to its AGN, NGC6951 also ex- 
hibits a pronounced SB ring at a radius of 5" (s480 pc) in Ha 
(Marquez & Moles 1993, Wozniak et al. 1995, Rozas et al. 1996; 
Gonzalez-Delgado & Perez 1997, Perez et al. 2000) and radio 
emission (Vila et al. 1990, Saikia et al. 1994 & 2002). Strong CO 
and HCN emission is associated with the SB ring (e.g., Kohno 
et al. 1999a, Garcfa-Burillo et al. 2005) while almost no emis- 
sion had been hitherto found in the centre of NGC 695 1 . Only 
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Table 1. HCN(l-O) line parameters, obtained at the various peak 
(i.e., not spatially integrated) by fitting a Gaussian profile to the 
(naturally weighted) data. Errors include uncertainties of the fit 
and calibration. " from Schinnerer et al., in prep. * from uni- 
formly weighted maps to avoid contamination by the ring emis- 
sion. c spatially integrated over the entire area of +9";"'inmJy; 
e in Jy km s _1 . 



recently have high angular resolution/high sensitivity PdBI ob- 
servations revealed faint CO(2-l) emission in the central 0.5" 
(Garcfa-Burillo et al. 2005; Schinnerer et al., in prep.). The latter 
observations are part of the PdBI NU(clei of)GA(laxies) project 
(e.g., Garcfa-Burillo et al. 2003). We observed NGC 6951 in 
HCN(l-O) to search for nuclear emission and assess differences 
between the SB ring and the AGN; the results of these observa- 
tions are presented here. 

2. Observations 

NGC 6951 has been observed at 3mm and 1mm using all six 
antennae of the IRAM PdBI in the new A, B, C and D config- 
urations (see http://www.iram.fr for telescope parameters) dur- 
ing January, February (A+B), April and May (C+D) 2006 . 
The phase reference centre was set to cnooo =20 h 37 m 14.470 s 
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Fig. 2. Iso- velocity map of HCN(l-O) (solid red & dashed 
blue lines) overlaid to the integrated HCN(l-O) emission {grey 
scale & grey contours (same as in Fig. [T}). Contours run from 
Vlsr=1349 km s _1 (dashed blue) to 1479 km s _1 (solid red) 
in steps of 10 km s _1 around the dynamical center (~10 km s _I 
off from the systemic velocity). The black lines indicate the PA 
(=160+20°) of the velocity gradient. 

and 5j2ooo =66°06' 19.70". The 3mm receivers were tuned to 
the frequency of HCN(l-O) shifted to the LSR velocity of 
vlsr =1424kms~ 1 , while the 1mm receivers were set to the 
13 CO(2-l) line. The 1mm data will be discussed in a separate 
paper (Krips et al., in prep.). Upper and lower sidebands at each 
frequency had 580 MHz bandwidth and 1.25 MHz resolution. 
Weather conditions were good throughout the observations with 
a water column varying between 4mm and 10mm and SSB sys- 
tem temperatures of 100-150 K (A+B) and 100-200 K (C+D) at 
3mm. 3C273, 3C454.3, 2200+420 and 1749+096 were used as 
bandpass calibrators while 1928+738 and 2037+511 were used 
as gain calibrators. Fluxes have been calibrated with MWC349 
and checked on 3C273, 1928+738 and 2037+511 with the flux 
monitoring program of the IRAM PdBI, resulting in an accuracy 
of ~10% at 3mm. The data were calibrated, mapped and ana- 
lyzed using the standard IRAM GILDAS programs CLIC and 
MAPPING. Using uniform weighting, the synthesized beam size 
is 1.25"x0.87" at a position angle (PA) of 72° at 3mm; natural 
weighting results in 2.78"x 2.27" at a PA of 1 12°. We reach an 
rms noise of ~0.8 mJybeairT 1 (~1.2 mJybearrT 1 ) at a spectral 
resolution of 5 MHz (=17kms _1 ) at 3.4 mm for natural (uni- 
form) weighting. 

3. Results & Discussion 

3. 1 . Starburst ring 

HCN emission is clearly detected along the starburst ring, con- 
sistent with Kohno et al. (1999a). The total flux observed with 
the PdBI is lower by a factor of ~2-3 than the flux obtained with 
the IRAM 30m telescope (Krips et al., in prep.), i.e., some of 
the HCN emission is resolved out by the PdBI. The two main 
maxima, which are seen in HCN by Kohno et al. (1999a) at an 
angular resolution of ~4.5" and peak at different positions than 



CO, split up into several sub-maxima (N,W,E,S; Fig.Q~|l at the 
higher angular resolution (<3") of the PdBI. The HCN peaks N 
and S are consistent in position with those in CO (Fig.Q~l Garcfa- 
Burillo et al. 2005). Additional HCN peaks are found to the west 
(W) and east (E) (see Fig.[TJ. The W-to-N and E-to-S peak ratios 
are higher in HCN (-0.6-1.0) than in CO (-0.4-0.6; Table Q} 
explaining why the positions of the merged peaks in Kohno et 
al. (1999a) do not agree between HCN and CO. The HCN-to- 
CO ratios (=Rhcn/co=Ihcn(i-0)/Ico(i-0); with 1= velocity inte- 
grated intensity) thus also differ between the S/N and W/E peaks. 
While Rhcn/co is found to be roughly ~0.03 at S and N (assum- 
ing Ico(2-i)/Icoa-0)=l; Table [U, we estimate Rhcn/co - 0.05 
at E and W. This might indicate a variation of the gas den- 
sity/temperature along the ring. The HCN kinematics in the ring 
agree well with those seen in CO and are hence not further dis- 
cussed in this letter. 



3.2. Central emission 

HCN(l-O) emission is detected in the central ~1" (^100 pc; 
=C) because of the higher sensitivity and angular resolution 
of our data compared to the one obtained by Kohno et al. 
(1999a), confirming previous hints found in CO(2-l) (Garcfa- 
Burillo et al. 2005; Schinnerer et al., in prep.). The nuclear 
HCN component C appears to be compact and unresolved in 
the PdBI beam. Assuming the "HCN conversion" factor of 
Xhcn(=Mhcn(H 2 )/Lhcn)=20^o M q (K kms" 1 pc 2 )" 1 from 
Solomon et al. (1992), we find a central dense gas mass of 
M gas (=M(H 2 +He)M2-lO)-lO 7 M which is a factor of -3-17 
higher than the mass of ~6-lO 6 M derived from the CO(2- 
1) line (assuming Ico(2-i)/Ico(i-0)=l (Table [T] see also Garcfa- 
Burillo et al. 2005). This large discrepancy cannot be solely due 
to the large uncertainties of the conversion factors but might 
further indicate non-standard gas conditions in the AGN. The 
HCN-to-CO ratio of the velocity integrated line flux amounts 
to >0.4 (Table Q] Ico(2-i)/Ico(i-o)=l)- This is a factor of -4-8 
larger than in the starburst ring and the value of 0.09 for the nu- 
cleus published by Kohno et al. (1999a) indicating that the nu- 
clear HCN is significantly enhanced in NGC 695 1 . Thus, unlike 
Kohno et al. (1999a), we conclude that the nuclear gas proper- 
ties in NGC 695 1 are not that different from those in the Seyfert 
galaxies NGC 1068 or M51, both of which have similarly high 
central HCN-to-CO ratios. The molecular gas chemistry in the 
AGN influenced regions of NGC 1068 and NGC 5194 is dom- 
inated by X-ray radiation (NGC 1068: e.g., Usero et al. 2004; 
M51: e.g., Matsushita et al. 1998) suggesting a similar scenario 
in NGC 695 1 . Thus, the central dense gas mass derived via HCN 
has to be considered an upper limit. 

The central HCN emission shows the steepest velocity gra- 
dient at PA»;160° with a velocity range of +70 km s _1 around 
the dynamical center over a radius of ~0.5" (Fig. [2}, eventu- 
ally indicating a gas rotation in a circumnuclear disk or torus. 
Even if accounting for the uncertainties of the PA, estimated 
to be +20°, the central kinematic axis is significantly different 
from the major axis of the galaxy (PA=135°) and of the bar 
(PA=100°). This non-alignment might be caused by a different 
inclination of the central gas disk/torus than of the galaxy, or, 
alternatively, by non-circular velocities or a warp. However, the 
central velocity gradient allows us to give a crude estimate of 
the enclosed dynamical mass. Assuming a radius of ~50 pc and 
a velocity range of +70km s , we find M dyn =5-10 7 -(sin i)- 2 M@, 
with /=inclination of the disk. This is similar to the gas mass of 
0.6-8-10 7 A4 Q , suggesting that the disk is seen closer to face-on 
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than to edge-on. Assuming further the inclination of the galaxy 
(~40°; Garcfa-Burillo et al. 2005) as rough estimate for the one 
of the central gas disk, the enclosed dynamical mass increases 
to ~2-10 8 /V(©. This is of the order of the nuclear black hole 
mass of ~2-10 8 A1o> estimated via the stellar velocity dispersion 
of the bulge (o-j=232km s _1 from Ho & Ulvestad 2001; see also 
Gebhardt et al. 2000). However, within a radius of 50 pc, stars 
might still contribute a significant fraction to the dynamical mass 
as well pointing towards an even lower inclination of the nuclear 
disk than assumed. 

4. Conclusions 

We have presented high angular resolution/sensitivity observa- 
tions of the HCN(l-O) emission in NGC 695 1 in this letter which 
exploit the new ABCD configurations at the IRAM PdBI. The 
main results can be summarized as follows: 

1. The two main HCN peaks in the ring seen by Kohno et 
al. (1999a) split up into several sub-maxima which coincide 
with those found in CO(2-l) by Garcfa-Burillo et al. (2005) 
and Schinnerer et al. (in prep.). However, the ratios of peaks 
within the ring differ between HCN and CO, probably indi- 
cating different gas densities/temperatures along the ring. 

2. We clearly detect compact HCN(l-O) emission in the nu- 
cleus of NGC 695 1 tracing a dense molecular gas mass in 
the range of ~0.6-6xl0 7 A1o. A position velocity cut taken 
along PA=(160±20)°, coinciding neither with the galaxy nor 
the bar major axis, yields the steepest velocity gradient in 
the nucleus. The gradient and the compactness of this cen- 
tral HCN emission suggest that it might arise in a rotating 
circumnuclear gas disk or torus with a radius of <50 pc. 

3. In contrast to Kohno et al. (1999a), we find significantly 
different HCN-to-CO line ratios in the starburst ring and 
the nuclear emission. While the starburst ring shows typi- 
cal Rhcn/co=0.02-0.05, HCN is significantly enhanced in 
the centre with Rhcn/co >0.4. Either the molecular gas 
in the nucleus is denser and/or hotter than in the starburst 
ring, increasing Rhcn/co, or the gas chemistry in the nu- 
cleus of NGC 6951 is dominated by X-ray radiation from 
the AGN, producing a higher nuclear HCN abundance and 
thus a higher Rhcn/co similar to the center of NGC 1068. 
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